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Abstract. This paper summarizes the limits on propagation of ultra high energy particles in 
the Universe, set up by their interactions with cosmic background of photons and neutrinos. 
By taking into account cosmic evolution of these backgrounds and considering appropriate 
interactions we derive the mean free path for ultra high energy photons, protons and neutri¬ 
nos. For photons the relevant processes are the Breit-Wheeler process as well as the double 
pair production process. For protons the relevant reactions are the photopion production and 
the Bethe-Heitler process. We discuss the interplay between the energy loss length and mean 
free path for the Bethe-Heitler process. Neutrino opacity is determined by its scattering off 
the cosmic background neutrino. We compute for the first time the high energy neutrino 
horizon as a function of its energy. 
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1 Introduction 

Observation of ultra high energy (UHE) particles, such as photons, ions and neutrinos, pro¬ 
vides the crucial information on astrophysical systems as well as mechanisms of charged 
particle acceleration in these systems. Such information cannot be obtained from the study 
of low energy emission, although it is much easier to detect. 

Propagation of UHE particles on cosmological distances involves interaction with other 
particles, as well as with electromagnetic fields, in the case of charged particles [1]. One of the 
most important reservoir of photons is the cosmic microwave background (CMB). Interaction 
with CMB imposes strong limits on propagation of UHE photons, protons, and nuclei. Ex- 
tragalactic background light (EBL), being the accumulated radiation in the Universe due to 
stars and active galactic nuclei, represents additional background of photons [2], which limits 
propagation of high energy photons. Yet another important background is cosmic neutrino 
background (Cz^B), which places a tight limit on the propagation of UHE neutrinos. 

In this work, we review stringent limits on propagation of UHE particles, namely pho¬ 
tons, protons and neutrinos, in the Universe due to their interactions with cosmic background 
of photons and neutrinos. We pay the particular attention to accounting for the cosmic evo¬ 
lution of CMB and Ci^B fields, being important at high redshifts. We discuss relation with 
previous results, as well as implications of new results, obtained in this work. 

First in Sec. 2 we discuss most important processes, responsible for interaction of UHE 
particles with cosmic backgrounds, as well as the corresponding cross-sections. Most of these 
processes are discussed in detail in [3]. Then in Sec. 3 we introduce the method used to 
compute the mean free path of UHE particles, which takes into account cosmological redshift 
of particle energy, as well as temperature evolution of the CMB and Cz^B. In Sec. 4 the 
definition of the mean energy loss distance is given. We present and discuss results in Sec. 5. 
Conclusions follow. 
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2 Processes 


2.1 Processes involving photons 

UHE photons are likely produced in sources of UHE cosmic rays. The most important 
process, responsible for intergalactic absorption of high-energy 7 -rays is the Breit-Wheeler 
process [4] for the photon-photon pair production 

71 + 72 — >-e'^ + e~ ( 2 . 1 ) 


It was first discussed by Nikishov [5] back in 1961 and then, after the discovery of CMB, by 
Gould and Schreder [ 6 ]. 

Breit and Wheeler [4] studied collision process (2.1) of two photons with energies E 
and £ in the laboratory frame, producing electron and positron pair. They found the total 
cross-section 
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h is Planck’s constant, rUe is electron mass, c is the speed of light and a is the fine structure 
constant. The necessary kinematic condition in order for the process (2.1) to take place 
is that the energy of two colliding photons is larger than the energetic threshold 2 meC^, 
i.e., X > 1. Due to this kinematic condition the function (2.2) has a low energy cut-off at 
X = 1. The cross-section has a maximum at x ~ 2, with ~ where or is Thomson 

cross-section. At higher energies it decreases as 1/x. 

A simple estimate of the mean free path for the Breit-Wheeler absorption of high energy 
photon can be given as follows. Considering the actual CMB photon density ncMB — 4:11 
cm“^ and taking (Tt /4 for the cross-section of the interaction, the mean free path is Xbw = 
{(TTncMB/^)~^ — 4.8 kpc. One can refer to this distance as to a horizon, namely the 
maximal distance to the source for which the particle with the given energy can still be 
detected on Earth. However, owing to the energy dependence of the cross-section, and 
cosmic evolution of the CMB photon field the actual mean free path strongly depends on 
energy. The characteristic energy of UHE photons interacting with CMB, having temperature 
today To ~ 2.725 K, is given by Ebw = {nieC^)/kT^ ~ 1.11 PeV. At lower energies, in the 
TeV range, photons interact by the Breit-Wheeler process with the EBL [ 6 - 8 ]. Hence the 
observation of TeV radiation from distant {d > 100 Mpc) extragalactic objects provides 
important constraints on the EBL [9-11]. 

At much higher energies the double pair production process 
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dominates, see [12, 13]. In this high-energy limit it has nearly a constant cross-section, see 
e.g. [3] ^ 

Clearly, this process has a threshold with the sum of energies of photons which must exceed 
4meC^. It imposes a limit for UHE photons propagation X^pp = {(ydppncMB)~^ — 121 Mpc. 
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2.2 Processes involving protons 

Charged UHE particles, such as protons and nuclei, are assumed to originate from extra- 
galactic sources, which work as ’’cosmic accelerators” [1]. Such particles interact with the 
CMB photons as well. In fact, the famous Greisen-Zatsepin-Kuzmin (GZK) limit [14, 15] 
was established by considering that this UHE particle interacts with the CMB photons via 
the pion photoproducton process 


P + 1 — 

and lose its initial energy. Due to the fact that at this process the proton loses more than 
half of its energy [16], such interaction imposes a strong cut-off on energies of UHE cosmic 
rays. The cut-off energy is easy to estimate. Recall that the characteristic energy in the 
Breit-Wheeler process (2.1) is Ebw = {rneC^)"^/kT^. When the photopion process (2.6) is 
concerned, the electron mass is exchanged with the pion mass, and an additional factor 4 
comes from the reference frame transformation, giving Ep^ = /kT^ ~ 3 x I^^Ebw = 

3.33 X 10^® eV. More careful evaluation of the energy by comparing energy losses due to 
photopion (2.6) and photoproduction of pair (2.8) processes (see below) gives the value 
Ep^ = 5 X 10^® eV [17]. The cross-section of the photopion process in high energy limit is 
constant with the value [16] 

Up^ ~ 120/i6. (2-7) 

The mean free path due to this process for energies E > Ep^ is = {crp^ncMB)~^ — 6 
Mpc. 

Another process relevant for interaction of UHE particles with CMB is the photopro¬ 
duction of electron-positron pair on a nucleus, or Bethe-Heitler process [18]. In the case of 
proton, which is the only one considered in this work, this process is 
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It has the characteristic energy Ebh = mempC^/{IkTo) ~ 1.0 x 10^^ eV. This process has 
a threshold with photon energy in the proton rest frame £' > 2meC^. We use for its cross- 
section in the proton rest frame the expressions given in [19], namely near the threshold with 
2 < e' < 4 
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where e' = £' j{meC^) is photon energy in the proton rest frame and jy = (e' ~ 2) / (e' -|- 2). 
At higher energies e' > 4 the cross-section is 
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where 5 = log(2e'). Expression (2.10) is logarithmically increasing at high energies, so we can 
take a characteristic value obtained by Bethe and Heitler aBH — (28/9)a [{ah)/(me c)]^ in 
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order to estimate the mean free path of UHE protons, which gives Xbh = {c^BHncMB)~^ — 
437 kpc. 

It is important to note that unlike the Breit-Wheeler process, leading to annihilation of 
UHE photons, or the pion photoproducton, where single interaction alters the energy of the 
UHE proton, the single Bethe-Heitler interaction does not change the proton energy signifi¬ 
cantly. Therefore, unlike all previous processes, the mean free path Xbh does not correspond 
to a horizon. Another quantity is used for this purpose, namely the mean energy loss dis¬ 
tance, defined as Xbh ~ [dE/{Ecdt)]~^, where E is the proton energy, which corresponds to 
the distance on which the energy of the UHE proton is reduced by a factor e due to numerous 
interactions with background photons [16, 17, 20]. However, it should be emphasized that 
single Bethe-Heitler interaction deflects the UHE proton by a small angle. This effect is 
discussed in detail below. 


2.3 Processes involving neutrinos 

UHE neutrinos can be produced either in astrophysical sources, or in some exotic new physics 
scenarios [21]. Below we compute the horizon due to interaction of UHE neutrinos with cosmic 
neutrino background (Ci^B). Eollowing [22] we assume Ci/B neutrinos are in their mass states. 
The cross-section is composed of two parts [22]. The resonant neutrino annihilation occurs 
in the s-channel; 

+ ^ f + f, (2.11) 

where bar denotes antiparticle, / is a fermion. It has a typical Breit-Wigner shape and is 
given in the analytic form in [23] . We take the small momentum expansion of the cross-section 
given by eq. (23) in [23] as 
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where Gp = 1.16637 x 10“^ GeV“^ is the Eermi’s coupling constant, ^ = (T/M^)^, T = 2.495 
GeV is the width of resonance and Mz = 91.1876 GeV is the mass of Z'^ boson, my is 
neutrino mass, E is energy of UHE neutrino in laboratory frame. Clearly, the position of the 
resonance scales inversely proportional to the neutrino mass. Throughout this paper we use 
the reference value m^ = 0.08 eV/c^, corresponding to the recent cosmological bound from 
the Planck mission [24], which gives characteristic energy E^ = Ml(^/2my ~ 5.2 x 10^^ eV. 
The amplitude of the resonance does not depend on neutrino mass, and is given by 


^fimax ^ 2V2GfMz/T ~ 0.471//6. 


(2.13) 


The second contribution is the non-resonant cross-section, which is adopted here in the 


form 
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where ~ 8.3 x 

We assume that neutrino are non-relativistic even at sufficiently high redshift, which is 
a good approximation for z < 10^ for = 0.08 eV/c^. Effects of non-zero momentum on 
the neutrino annihilation cross-section are studied in [22, 23]. Using the number density of 
relic neutrinos ncuB — 112 cm“^ and the non-resonant cross-section in the high energy limit 
one can estimate the horizon for UHE neutrinos at highest energies. Solving the Friedmann 
equation (see next section) one finds for the redshift Zy ~ 84. 
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3 The optical depth and the mean free path 


In this section we compute the optical depth for the propagation of UHE particles in the 
Universe. Imposing the condition that it equals unity we determine the corresponding mean 
free path. It should be noted that simple estimates, made in the literature, as well as in 
previous section, do not take into account evolution of CMB and Cz^B fields with time. The 
simplest way to account for cosmic redshift is to compare this estimate of the mean free path 
to the expansion scale c/Hq, where Hq is the present day Hubble parameter, see e.g. [17, 20]. 
In what follows we describe more rigorous way to take into account both redshift of particle 
energy as well as the evolution of CMB and Cz^B fields with redshift. 

The optical depth along the particle world line C is defined as 

r = j^crj^dx^, (3.1) 

where cr is the cross-section of a given process, is the 4-current of particles, on which the 
UHE particle scatters, and dx^ is the element of the UHE particle world line. We assume 
the Universe is homogeneous and isotropic, and the background particles are either CMB 
photons or Cz^B neutrinos. Both have thermal distribution functions, given by 

( 3 - 2 ) 


where k is the Boltzmann constant, T is the CMB or Cz/B temperature, the sign ” is 
for photons while the sign is for neutrinos, £ and /i are the energy and the chemical 
potential of background particles (for photons fj, = 0). Then the optical depth (3.1) is 

riE,t) = l^dt' £^d£fi£)a{E,£,t'), (3.3) 

where Etr is threshold energy in a given process, = 2 is the number of helicity states 
for both protons and neutrinos. Here we assumed that UHE particles move along light-like 
geodesics. The integral over time can be transformed into the integral over redshift by means 
of the Eriedmann equation. The latter for the flat Universe reads 

1 da\^ 8 ttG . 
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where a is the scale factor, p is energy density of the Universe, G is Newton’s constant. Erom 
this equation, the definition of cosmological redshift, as well as the definition of the density 
parameters 

ao/a = l + z, £1.1 = —, = (3.5) 

Pc OTTCt 

where Hq and oq are present time Hubble parameter and scale factor, respectively, we have 
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where Hq is the Hubble parameter and the function H{z) is given by 
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and Qr, and r^A are present densities of radiation, matter and dark energy, respectively. 
Then the expression (3.3) can be written as follows 


t{E,z) = 
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Cosmic expansion results in the energy and temperature dependence on redshift 
T = (1 + z)Tq, £ = {1 + z)£q, E = {1 + z)Eq^ 


(3.8) 


(3.9) 


where temperature ~ 2.725 K for photons, = (4/11)^^^ ~ 1.95 K for neutrinos and 
energies Eq,£q are measured at the present time. 

The second integral in (3.8) can be simplihed, provided two conditions are fulfilled; a) 
the cross-section does not depend on the energy of background particle and b) there is no 
threshold in the given process {£tr = 0). In this case one has 

£'^d£ f{£)a{E, z) = criE, z)n(z) = cr{E, z)no (1 + z)^ , (3.10) 

where no is present number density and it stands for either 
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for photons, or no,i/ = 3/4:{To^y/To^^) ~ 113 cm for neutrinos. Then eq. (3.3) becomes 
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When the cross-section is just a constant, the integral (3.11) can be readily performed. 
Assuming — 9.2 x 10“^, £Im — 0.315, ~ 0.685 and Hq = 67.3 km/s/Mpc [24] in the 

matter dominated epoch we have 



{1 + z'fdz' 

[nMii + z'r + £iAV/^ 




1.045z, z <C 1, 

1.006z3/2, z > 1. 


The mean free path is defined by the condition t{Eq,z) = 1. For the constant cross- 
section a at low redshift 2 : ^ 1 we get the traditional dehnition A = {an)~^ used above. For 
high redshift z ^ 1 one can dehne the redshift, corresponding to the mean free path as 
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Using this equation we obtain for UHE neutrinos with highest energies za — 84. 


4 The mean energy loss distance 

When UHE particle annihilates in a given process, such as in the case of Breit-Wheeler one 
(2.1), the mean free path correpsonds to the horizon dehned above. 

Another possibility is that the particle is not annihilated in a given process, but scat¬ 
tered, such as in the case of proton producing the pion (2.6). When the energy loss in single 
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scattering corresponds to a large fraction of UHE particle energy, the situation is similar to 
the case of annihilation. However, UHE particle may lose only a small fraction of its energy, 
as in the case of Bethe-Heitler process (2.8). Here another relevant quantity corresponds the 
particle horizon dehned above is the mean energy loss distance A. We define it following [25] 
as 


A-i 


\ E cdt) 
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Then we evaluate the quantity 


T = 
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It is computed below for the Bethe-Heitler process. 


5 Results 


Now we apply the method developed in the previous section to the computation of the mean 
free path for UHE photons, protons and neutrinos, as well as the mean energy loss distance 
for protons interacting via the Bethe-Heitler process. 


5.1 Photons 


First, we consider cosmic limits on propagation of UHE photons. In the Breit-Wheeler 
process (2.1) the cross-section depends on both energies through the definition (2.3). When 
one considers all possible orientations of CMB photons additional averaging over their angular 
distribution has to be performed [5, 6]. The resulting averaged cross section differs from eq. 
(2.2). The useful approximations for this quanity can be found e.g. in [6, 26, 27]. We use 
the accurate expression given by eq. (3.23) in Ref. [1]: 

(:, + il„gI_i + T)x (5.1) 
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We change our variables in eq. (3.8) and get 
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and yo is the energy Eq in units of the critical energy Ebw = /kTo ~ 1.11 x 10^® 

eV. The intergal over energy can be evaluated numerically and we find a reasonable fit 


Fi (y) = 0.839 (y^'^ + 2x 10“®y2.8^ 



(5.5) 
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Then eq. (5.2) becomes 


A 1 rlr' 

r,,(!/0, ^ / (inr flT) ^ 

This integral is also evaluated numerically. 

In the low energy E <C Ebw and high redshift z 3> 1 limit the integral (5.6) can be 
evaluated analytically. The redshift corresponding to the mean free path in this limit is 


z\BW - 0.21 



(5.7) 


This result is known as the Fazio-Stecker relation [28], see their eq. (9). 

In addition to the Breit-Wheeler process (2.1), following [13] we consider also the double 
pair production process (2.4) with the cross-section defined in (2.5). This process is relevant 
for the highest energies. The optical depth for this process is 
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(5.9) 


and B Ri 15.3. 

The condition T{yo,z) = T^^{yo,z) + Tyy^{yo, z) = 1 in eqs. (5.2) and (5.8) determines 
the mean free path of UHE photons. This mean free path is shown in Fig. 1 in megaparsecs 
and in Fig. 2 in cosmological redshift by the solid curve. The region above the solid curve is 
opaque for high energy photons. In addition, thick black dashed line shows the boundary of 
transparency for EBL, according to the baseline model of Inoue et al. [29], while dotted line 
in Eig. 1 shows the cosmological horizon z = 00 . 

For the distance smaller than a critical value of about dc = 6.8 kpc, the CMB is 
transparent to high-energy photons with arbitrary energy. For larger distances there are 
two branches of solutions for the condition T^^{yQ,z) = 1, respectively corresponding to 
the different energy-dependence of the average cross-section (5.1). This average cross-section 
d'.y.y (x) increases with the center of mass energy x from the energy threshold x = 1 to x ~ 3.5, 
and decreases from x ~ 3.5 to x — 00 . The energy of the UHE photon corresponding to 
the critical distance dc is about 1.11 PeV, which separates two branches of the solution. The 
double pair production process (2.4) is relevant for the highest energies, as expected. Photons 
with energies above 10 PeV are absorbed by the double pair production if they are emitted 
at redshift above z ~ 0.03 (distance about 120 Mpc). 

For comparison we show by the dashed curve also the mean free path computed at 
z = 0, namely neglecting cosmological expansion, see e.g. [13]. 


- 8 - 













d, Mpc 



E, eV 


Figure 1. The mean free path, measured in megaparsecs as a function of energy E of UHE particles, 
measured in electronvolts. In the region above the curves the optical depth is larger than unity. 
Thin black curve shows the mean free path of UHE photons. Dashed black curve shows the photon 
mean free path computed without accounting for cosmological evolution (imposing z = 0). Thick 
blue dashed curve shows the boundary of transparency for extragalactic background light (EBL), 
according to the baseline model of Inoue et al. [29]. Thick blue curve shows the mean free path 
of UHE protons (GZK limit). Blue dashed (dotted-dashed) curve shows the mean free path (mean 
energy loss distance) for UHE protons due to Bethe-Heitler process. Dotted horizontal line shows 
cosmological horizon. 


Z 



E, eV 


Figure 2. The same as in Fig. 1 for the distance measured in cosmological redshift. Red thick curve 
shows the mean free path of UHE neutrinos. 


We also show by the dashdotted thick curve the mean free path for the photon-photon 
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scattering which follows from the Euler-Heisenberg lagrangian, see e.g. [3, 30]. We will 
discuss this process in a separate publication [31]. 

Finally, the black dotted thick curve shows the horizon of photons with energies above 
20 GeV and below 100 TeV, which is determined by their interaction with the EBL. The 
latest EBL model [29] is used. It is clear, that the contribution of CMB photons gives the 
absolute upper limit on the mean free path. In the energy range between 1 GeV and 20 GeV 
the propagation of high energy photons is limited only by the CMB radiation. 


5.2 Protons 

Second, we consider the propagation of UHE protons, accelerated in a source located at a 
cosmological distance from Earth. First, considering the photopion process (2.6) we use the 
method developed in the previous section and compute the GZK limit [14, 15]. This limit 
applies to protons and other charged particles, leading to the existence of a cutoff in the 
observed spectrum of (UHE) cosmic rays at about 10^® eV. For the photopion process (2.6) 
one can use the simple expression (3.11) with the constant cross-section (2.7). However, we 
compute the optical depth in the same way as in the case of double pair production, using 
eq. (5.8) with different value of the constant B' ~ 253. 

The mean free path due to photopion process is shown by the blue thick curve in Fig. 
1 in megaparsecs and in Fig. 2 in cosmological redshift. From Fig. 1 it appears that for 
energies well below Ep^ ~ 3.3 x 10^° eV the GZK limit approaches the cosmological horizon. 
Instead, from Fig. 2 it follows that the mean free path measured in redshift, below the energy 
increases with decreasing energy as a power law, which is a consequence of eq. (3.12). 
Similarly to the Breit-Wheeler case, the integral (5.8) is evaluated analytically in the 
low energy E <C Ep^ and high redshift z S> 1 limit, with the result 


zx.GZK - 0.57 




(5.10) 


We also evaluate the mean free path due to the Bethe-Heitler process (2.8). Since cross- 
sections (2.9) and (2.10) are given in the proton rest frame, one has to transform photon 
energy to this reference frame using 


T^C 

£’ = 2T£ =2-(5.11) 

mpC^ 

where the primed quantity corresponds to the proton rest frame, while unprimed quantities 
to laboratory reference frame. Then it is convenient to make use of the same type of variable 
change as before for the Breit-Wheeler process^, with a difference that instead of electron 
mass squared a product of electron and proton masses arises, namely 


X = 2 


E£ 

meUipd^ ’ 


- _ r, ^0 kTo _ Eq 

^ 2 2 TP 

nipC^ rUeC^ Ebh 


(5.12) 


^We assume that UHE protons collide with the with CMB photons head on. More accurate calculation 
with average over angular distribution of the CMB photons does not change qualitatively our results. 
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The optical depth is computed in the laboratory frame as follows^ 
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The intergal over energy can be evaluated numerically and we find a reasonable fit 
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Then eq. (5.2) becomes 
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(5.16) 


The mean free path for protons interacting via the Bethe-Heitler process is shown by blue 
dashed curve in Fig. 1 in megaparsecs and in Fig. 2 in cosmological redshift. The integral 
(5.15) is evaluated analytically in the low energy E <C Ebh and high redshift z S> 1 limit, 
with the result 

/ ^ \-i /2 

. (5.17) 


Z\,BH ^ 0.43 


\Ebh) 


As discussed before, at energies above 10^® eV the mean free path is relatively small, 
about a few hundred kiloparsecs, and it quickly decreases with increasing energy. This is in 
contrast with the large mean energy loss path, which is above 1 Gpc at energies 10^® — 10^*^ 
eV, see e.g. [17, 20]. It means, that before UHE proton starts to lose its energy, it is scattered 
many hundred times [16]. At each of this scattering the proton recoils, being deflected by a 
small angle measured in the laboratory reference frame. From the analysis of the cross-section 
as the function of recoil angle [32], see also [33] in the proton rest frame, it follows that in 
the high energy limit the photon recoils in the plane, orthogonal to the incident photon. It 
implies that each scattering produces a deflection of the UHE proton in the laboratory frame 
by angle ~ I/ 7 , where 7 = E/ (mpC^) is proton Lorentz factor. The number of scatterings 
is given approximately by r. 

One can compute the mean energy loss distance defined in (4.1) and then evaluate the 
quantity (4.2) using eq. (18)-(20) in [25] to obtain 


T = 


f 


dz' 


D 


y^Jo {l + z')^H{z') 


-f 

') J2 


dx 


where the function (j) (x) is given in eq. (16) in [25] 


exp(f) - 1 


(p{x) = x -86.07 + 50.95 log x - 14.45 (log xf + 2.667 (log x) 


(f){x) 


;^\3 


(5.18) 


(5.19) 


^Note that the integral over energy is not transformed to the proton rest frame, as done e.g. in [25]. 
Instead, only a change of variables is performed in the integral (5.13). 
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Figure 3. The optical depth at the mean energy distance for the Bethe-Heitler process. 


and 


TT^ nipHo \mcj \meC^J 


(5.20) 


From the condition f = 1 we determine the mean energy loss distance A. This distance is 
shown by blue dash-dotted curve in Fig. 1 in megaparsecs and in Fig. 2 in cosmological 
redshift. 

We evaluate the optical depth (5.13) at the redshift, corresponding to A for energies in 
the range between 10^® and 10^^ eV and find it in the range 10^ to 10^, see Fig. 3. Since the 
deflection at each interaction is small, the average number of interactions is proportional to 
the optical depth. The average deflection angle is then 


7 


(5.21) 


We hnd an average deflection angle of UHE protons as a function of proton energy for sources 
located at the mean energy loss distance and show it in Fig. 4. At energies E = 10^®eV, this 
angle is about 5 ~ 15'' and it decreases down to 5 ~ 2.4 mas for E = lO^^eV. This latter 
scale is larger than the angular size of distant blazars [34]. Such deflection, although small 
compared to deflection in galactic magnetic field [35], inevitably leads to dimming of point 
sources of UHE protons. It makes also more difficult to detect distant sources. 


5.3 Neutrinos 

Third, we consider the propagation of UHE neutrinos. Such neutrinos can be produced in 
the source of UHE cosmic rays in decay of secondary pions tt'*' —)• or secondary 

neutrons n —?> p -|- e“ -|- Fg [16]. UHE neutrino can be produced also in some extensions of 
the standard model of particle physics [21]. Such UHE neutrino interacts with the Ci^B via 
the process (2.11). 
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Figure 4. The average deflection angle of UHE protons as function of proton energy for sources 
located at the mean energy loss distance for the Bethe-Heitler process. 


The cross-section of this process has a resonance, and it approaches a constant for 
highest energies. We compute the optical depth, which instead of eq. (3.8) is given by 


Ty^{E,z) = 


1 


Hq 


f 


dz' 


nWHoJo {l + z')H{z')Je,^ 
1 c 

-no,u I 

0 


— {myC^)^d£f{£)a{E,z) « (5.22) 


{l + z'fdz' 
H{z’) 


a{E{l + zi)), 


using the cross-sections given in the laboratory reference frame by eqs. (2.12) and (2.14). The 
mean free path for neutrinos, measured in cosmological redshift, is shown in Fig. 2 by the 
thick red curve. Since the characteristic redshifts are high, this curve practically coincides 
with the horizon, when measured in megaparsecs, so we do not show it in Fig. 1. It is 
clear that the Breit-Wigner resonance in the cross-section decreases the mean free path in a 
wide range of energies. The lowest redshift for E zz Ej. at which the Universe is transparent 
for UHE neutrinos is Zmin — 30. The resonance produces a dip around Er/{1 + Zmin) — 
1.7 X 10^^ (m,y/0.08 eV)~^ eV, where Zmin — 30. Additional broadening of the resonance, 
due to thermal effect, is discussed in detail in [22]. At higher energies the corresponding 
redshift is z ~ 87. 

Similarly to the previous cases, in the low energy E <C Ep^ and high redshift z ^ 1 
limit, we find 

/ E\ 

- 14 j . (5.23) 


6 Conclusions 

We reviewed cosmic limits on propagation of ultra high energy particles such as photons, 
protons and neutrinos, set up by their interactions with the cosmic background of photons 
and neutrinos. In doing so we take into account explicitly cosmic evolution of both cosmic 
backgrounds, and redshift of UHE particle energy. This is in contrast with majority of the 
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literature, where corresponding mean free paths are found at present epoch, neglecting cosmic 
expansion. A number of new results were obtained, in particular: 


• for UHE photons the contribution of CMB photons gives the absolute upper limit on 
the mean free path. At high redshift, where other radiation backgrounds, such as EBL 
are absent, the CMB radiation limits the propagation of UHE photons at energies above 
GeV. 

• for UHE protons the mean free path due to Bethe-Heitler process appears to be much 
shorter than the mean energy loss distance. This results in multiple deflections suffered 
by UHE protons, before they start to lose energy in the energy range 10^® — 10^® eV. 
Such deflections result in dimming of point sources of UHE protons, which makes it 
more difficult to detect them. 


• for UHE neutrinos for the first time we compute the horizon as a function of redshift. 
We found that the Universe is transparent of UHE neutrinos at redshifts z < 30, near 
the Breit-Wigner resonance at Er — 5.2 x 10^^ {rriu/O.OS eV)~^ eV, and it is transparent 
at redshifts z < 87 at higher energies. 


Remarkably, in the low energy and high redshift limit, the Eazio-Stecker relation [28] 
holds for all processes with exception of neutrinos, and it is given by a universal ex- 

pression — O (1) ; where is the characteristic (e.g. threshold) energy 

for a given process. In the case of neutrinos similar power law exists za oc E~‘^I^ . 
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